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SUMMARY 



The existence of energy dissipations in gas dynamics, 
which must be attributed to a lag in the vibrational heat 
ca;oacity of the gas, has been established both theoreti- 
cally and experimentally. 

The flow about a very small impact tube is discussed. 
It is shown that total-head defects due to neat-capacity 
lag during and after the compression of the gas at the 
nose of an impact tube are to be anticipated. Experi- 
ments quantitatively verifying these anticipations in 
carbon dioxide are discussed. A general theory of the 
dissipations in a nope general flow problem is developed 
and applied to some special cases. It is pointed out 
that energy dissipations due to this effect are to be 
anticipated in turbines. Dissipations of this kind might 
also introduce errors in cases in which the flow of one 
gas is used in an attempt to simulate the flow of another 
gas. Unfortunately, the relaxation times of most of the 
gases of engineering importance have not been studied. 

A new method of measuring the relaxation time of 
gases is introduced in which the total-head defects ob- 
served with a specially shaped impact tube are compared 
with theoretical considerations • A parameter is thus 
evaluated in which the only unknov.ni quantity is the re- 
laxation time of the gas. This method has been applied 
to carbon dioxide and has given consistent results fcr 
two impact tubes at a variety of gas velocities. 



INTRODUCTION 



The heat content of gases is primarily three forms 
of molecular mechanical energy. .First, there is the 

translat ional kinetic energy which is TjrRT, where R 
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is the gas constant and T is the absolute temperature. 
Secondly, there is the rotational kinetic energy. For 
all gases near or above room temperature, the n rota- 
tional degrees of freedom involving moments of inertia 
due to the separation of atomic neuclei have energy states 
close enough together that the rotational internal energy 

is close to the classical value ^RT. The third prin- 
cipal form of internal energy is the vibrational energy 
of the molecule s 0 If the frequencies of the normal 
modes of vibration of the molecule are known (say, from 
spectra), the vibrational heat capacity can be computed 
by the methods of statistical mechanics. (See, for 
example, reference 1.) 

The possibility of dispersion and absorption of 
sound due to parts of the heat capacity lagging behind 
the rar)id temperature changes accompanying the propaga- 
tion of a sound wave In a gas fas first discussed theo- 
retically by Jeans and Einstein • Dispersion and ab- 
sorption In carbon dioxide observed by Pierce were shown 
by Kerzfeld and Rice to be attributable to lagging of 
the vibrational heat capacity of the gas. Kneser was 
able to account quantitatively for dispersion and 
absorption in COo and oxygen on the assumption that the 
vibrational heat capacity lagged. 

The dispersion and absorption of sound in several 
gases have been investigated and a fairly complete 
bibliography is available in reference 2. It is found, 
in general* that dispersion and absorption many times 
larger than those attributable to vicosity and heat 
conduction are to be expected in gases With vibrational 
heat capacity. These effects can be describee? by rela- 
tions such as those given by Kneser and can be attributed 
to the vibrational heat capacity of the gas. 

All the measurements of dispersion and absorption 
have demonstrated that most impurities markedly reduce 
the relaxation time of a gas; for example, Kneser and 
Enudscn (references 3 and 4) concluded that the adjust- 
ment of the vibrational heat capacity of oxygen was 
dependent entirely on the action of Impurities. 
Various experiments with COo have shown that, at room 
temperature, collisions with water molecules are 500 | 
times as effective as collisions with C0 2 -molecules in 
exciting vibration in COg-molecules . This strong 



dependence on purity has produced great discrepancies 
among the relaxation times measured by the various workers 
in this field. There has been inuch better agreement 
among the measurements of the effectiveness of impurities. 
Sonic measurements in CO^ are discussed in appendix iu 
a quantum-mechanical theory of relaxation times developed 
by Landau and Teller is discussed in appendix 3. Two 
conclusions, which are verified by the sonic work in CO2, 
are important to the present paper: (1) All the vibra- 
tional states of a single normal mode adjust with the same 
relaxation time and (2) the logarithm of the relaxation 
time (expressed in molecular collisions) is proportional 

to T °. Verification of conclusion (2) is presented 
in figure 1. 

Dr. Vannevar Bush helped to initiate this work by 
asking the writer a stimulating question. The author 
also is very grateful to Professor Teller for helpful 
discus s ions . 



EFFECTS IN GA3 DYNAMICS 



In the flow of gases about obstacles, compressions 
and rarefactions accompanied by temperature changes 
occur* The time in which these temperature changes 
take place is controlled by the dimensions of the ob- 
stacles and the velocity of flow. If these time inter- 
vals are comparable with or shorter than the time re- 
quired for the gas to absorb its full heat capacity, 
the gas will depart from its equilibrium partition of 
energy. In this case, the transfer of energy from parts 
of the heat capacity that have more than their share to 
parts that have less than their share will be an irre- 
versible process and. will increase the entropy of the gas. 
If the time intervals involved are comparable with the 
relaxation time of the gas, this increase in entropy can 
be used to measure the relaxation time of the gas (refer- 
ence 5) . 

Turbine -working fluids such as steam, air, and ex- 
haust gas have appreciable vibrational her.t capacity at 
high temperatures, If these gases have relaxation times 
comparable with or shorter than the intervals during 
which temperature changes occur in the gas, losses 
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attributable to heat-capacity lag must be anticipated. 
A. rough estimate has indicated t&kt, unless the relaxa- 
tion times of the working fluids are very short, the 
losses at high temperatures due to heat-capacity lag can 
be comparable with the losses due to skin friction. 
Unfortunately | no measurements .of the relaxation times 
of the usual turbine -working fluids exist. 

Various persons have proposed, in wind-tunnel tests 
and in tests of rotating machinery $ the substitution of 
gases that have properties enabling tests to be made more 
conveniently at a given Hach number or Reynolds number 
than wi 4 :}.! the actual working fluid. In such cases, 
care must be taken to ensure that an error due to dif- 
ferences in heat-capacity-lag behavior of the fluid used 
and the working fluid is' not introduced. For example, 
according to a rough calculation, a. wing in pure CO 2 
might have a drag coefficient twice as large as the same 
winr in air at the sane Ilach number and Reynolds number. 

In the fallowing discussion, the existence of these 
dissipations in gas dynamics is lemons trated and a gas- 
dynamics method of measuring the relaxation times is 
developed. The application of this method to the meas- 
urement of the relaxation tines of gases of engineering 
importance is proposed. 

FLOV; ABOUT A VERY SMALL IMFACT TUBS 



As a first example of the energy dissipations to be 
expected from heat-capacity lag, consider the total head 
measured by an impact tube' in a perfect gas. For defi- 
niteness, consider the apparatus illustrated schematically 
in figure 2. The gas enters the chamber and settles at 
the pressure p 0 and the temperature Tq. It then ex- 
pands to a pressure pn and a temperature Tj out of the 
faired orifice, which is designed to give a temperature 
drop gradual enough that the expansion through the orifice 
is isentropic. The gas that flows along the axial stream 
line of the impact tube is then brought to rest at the 
nose of the tube and, during this process, its pressure 
rises to p 2 and its temperature rises to To. If this 
second process is slow enough to be isentropic also, the 
entropy and the energy of the gas that has reached equi- 
librium at the nose of the impact tube are equal to the 



5 



values in the chamber and hence the pressure Pg equals 
p 0 and the reading on the alcohol nanometer Is zero. 

Con.33.der, however, the other extreme case in which 
^ tho compression time - that is. tho time required for the 

-sf gas to undergo the greater part of its temperature rise - 

at the no se of the impact tube is small compared with the 
time required for the gas to absorb its full heat capacity. 
The orifice la considered largo enough that, during the 
expansion through it, tho gas maintains equilibrium. A 
pkrt of the heat capacity of the gas does not follow 

the rise in temperature during the compression as the 
gas is brought to rest at tho nose of the impact tube and 
adjusts irreversibly after the compression is over. The 
resultant increase of entropy in this case means that the 
pressure p ? is lower than p 0 . This increase in entropy 
is now calculated© All temperature changes are assumed 
small enough that the heat capacities of the gas can be 
taken as constants 

At the beginning of the adjustment, the lagging part 
o^ the heat capacity e^-u is still in equilbrium with 
a thermometer at the temperature^ fr, -| while the transla- 
tional and other degrees of freedom with heat capacities 
totaling Cr> ! , the relaxation time of which can be neg- 
lected, are" in equilibrium with a thermometer at some 
higher temperature T. Energy then flows from the heat 
capacity c p ! to the" heat capacity, c^-^, increasing 

the temperature TrHt) associated with c v -b from T l 
to the final equilibrium, temperature, v/hich is T Q • 
Conservation of energy gives the following relation 
between T and -yi^ * • 

c vib T vib + V T = C P T 0 (1) 



where c p is the total heat capacity at constant pres- 
sure. • . The entropy increase when an. element of energy 
dq- flows from T to Tyife Is 

j q _ dq dq / 1 l\ [n) 
dS - « ~ c Tr *>. d? ... — ) 

VxD V10 1 rp rp J 



i vib . - \ x vib 
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Eliminating T in equation (2) from equation (1) and in- 
tegrating over the whole process gives 



A3 = 



/ -V + i f 

j c Tlb . c vtt T vlb Erf*, 




- leg "!»» -1 U2 (3) 



Equation (3) gives the entropy difference between 
the gas in the chamber and the gas at equilibrium at the 
nose 'of the impact tube. Because the energy and hence 
the temperature is the same at the beginning and at the 
end of the process, the ratio of chamber pressure to 
impact-tube^ pressure can readily be computed from the 
perfect gas relation 

S r c p log T - R log p + Constant 

which gives 

PO 

AS - R log — 

and 



PO 
P2 




(4) 



It nay be instructive to derive this relation by 
considering the isentropic parts of the process. During 
the slow expansion, the enthalpy theorem (see appendix C) 

O 

gives c^T + iu^ = Constant (u is flow velocity) and, 
during the instantaneous compression, Cp 1 ! + tju = Constant, 
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°vib teing omitted because it taker, no part in the com- 
pression. Combining these equations gives 



LT\ 



(5) 



where Tg is the temperature reached by the translational 
decrees of freedom before the adjustment period starts. 
The adlabatic-compression relation can be applied to both 
the expansion and. the compression wltja the appropriate 
heat capacities to calculate pgj thus, 



To 



r Pc \ c p 




T l 



(6) 



Combining equations (5) and (6) gives, after several 
manipulations , 




vib 



which is of course identical with equation (4)^. 

Pn - P2 

The percentage total-head defect ICO is 

" ' , ?0 - Pi 

plotted against chamber pressure Pq/Pi * n fi£U" e 3 for 

Cp ? = 3.5R. The apparatus schematised in figure 2 was 
used to check equation (4) for COo where the vibra- 
tional heat capacity would be expected to Xag« The 

orifice was a hole in a 75- inch plate with Its diameter 
variation designed for constant time rate of temperature 
drop. The last l/l6 inch of the flow passage was 
straight in orde^ that the streamlines in the jet would 
be straight and axial and hence the static pressure at 
the orifice exit would equal atmospheric pressure. The 
glass impact tube was 0.005 inch in diameter and its end 
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were between 300 and 600 feet nor second. The expansion 
therefore took place in times ranging between 1*4 x 10 x 
and 2.8 x 10""* second* The compression at the nose of 
an impact tube takes place while the gas flows a distance 
of the order of 1 tube radius. (See fig. 5*) The com- 
pression times then ranged between 7 X lO" 7 and 
14 x 10*' second. Commercial CO2 was used and, be- 
cause it was fairly dry, a relaxation time of the order 

of l.C~ 5 second was expected. It seemed likely, therefore, 
that this setup would approach the case of an r'sentropic 
expansion and an instantaneous compression closely enough 
for the results to bear at least a qualitative resemblance 
to equation ( 4) . 

Preliminary to the investigation of heat-capacity 
lags, it was necessary to make sure that hydrodynamic 
effects other than heat-capacity lag would not produce a 
reading on the alcohol manometer. Air and later nitrogen 
at room temperature were therefore substituted for COg 
at the beginning of each run. It was always found in 
these preliminary tests that, v.hen the tube was properly 
alined, the difference in pressure measured by the alcohol 
manometer was very small and could be accounted for 
entirely by lags in the small vibrational heat capacity 
of sir (about 0.02R) . 

Carbon dioxide was then introduced into the apparatus 
and the observations shown in figure 4 were made. The 
gas was heated before entering the chamber, and its 
temperature was measured by a 3mall ' thermocouple inserted 
in the jet close to the impact tube. In accordance with 
aerodynamic experience, the temperature measured by the 
thermocouple was assumed to be 0.9T 0 + 0.11V The dlf ~ 
ference between Tq and was always less than 50° P, 

corresponding to a difference in c vi - b of less than 
8. percent, and was thus considered accurate enough to 
assume a constant c v i b and to compute this value at a 

m 4-"rp 

temperature T E — J* # 



The pressure p Q - p« was read by the mercury 
manometer, p-j^ by a barometer, and p Q - p 2 by the 
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alcohol nanometer, which was fitted with a microscope to 

make possible readings to 0.001 inch. 

In figure 4 the reading of the alcohol manometer is 
plotted against the chamber pressure Pq/pi* The ex- 
perimental values at both temperatures agree with the 
^ theoretical values more closely than could have been 

anticipated* It will become clear later that the theo- 
retical and experimental values agreed so closely becar.se 
small entropy increases in the orifice, attributable to 
too-rapid expansion, just about compensated for the fact 
that the compression was not quite instantaneous compared 
With the relaxation time of the gas. It should be 
pointed cut that ordinary hydrodynamic effects such as 
misalinement of the impact tube would be expected to pro- 
duce a total-head defect which would vary directly as 



Po 



Pi 



- 1 



GENERAL THEORY OP ENERGY DISSIPATIONS IN GASES 
EXHIBITING HEAT -CAPACITY LAG 



In the general case in which the temperature changes 
may be neither very fast nor very slow compared with the 
relaxation time of the gas, the temperature history of a 
gas particle as it flows along a streamline must be con- 
sidered. The problem is greatly simplified if the effect 
of heat-capacity lag on velocity distribution is 
neglected in order to get the effect of the lag on energy 
dissipation* This procedure can be regarded as the 
first step in an iteration process and is probably ade- 
quate for the applications now contemplated. The re- 
striction that the temperature changes involved in the 
flow arc small enough for the heat capacities to be con- 
sidered constant is also retained* 

Assume, therefore, that the velocity distribution in 

the field of flov; is determined by standard gas-dynamics 

methods* The velocity distribution is usually given as 

a function of space coordinates u(x,y*z) or along the 

^-streamline as u Q (s) , where s is the distance along 

the streamliner This expression can be converted to a 

function of time u -(t) by integration of dt = J&£ — 

° u(s) 
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along a streamline • The function u a (t) is taken for 
granted and the entropy increase in the flow along a 
streanline is determined* 

By introducing the variable c, which represents 
the excess energy per unit mass in the lagging heat 
capacity over the energy at equilibrium partition at the 
transnational temperature T, it is seen that 

c D T + -u 2 + e~ Constant (7) 

The assumption is now introduced that there is only 
one type of heat energy in the gas E 4^ which lags 
appreciably behind the translation temperature and that 
its time rate of adjustment is proportional to its de- 
parture from equilibrium; that is, 

dt 

This assumption is in agreement with the spnic theories 
previously discussed. From the definition of € , 



* - e - c , : 

vib vib 



because c v ik T v -?"b is ^ e ^QUiHb^um value of E vib 
(measured from an arbitrary zero) • By combining these 
equations, s v ib can " oe e i iril i nated - to yield 



de 



dt V 



- -c_. v . % ~ k€ (8) 



Ciu 



The meaning of k can be made clear if the variation of 
c with time is examined for the case in which the total 
heat energy of the gas remains constant. In this case, 

c T + e » Constant 
-o 
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Equation (8) then becomes 



<?€ _ c vib d€ 



dt 



or 



dt " Crj f 

Cr) 

from which k— ~ is the reciprocal of the relaxation 
C P 

time t of the fas. It will te seen that these equa- 
tions are restricted to gases with only one relaxation 
time • 

In order to simplify later expressions and to 
clarify their physical meaning, there are introduced the 
dimensionless variables 



-vib 



K a 



t 

I7/U 



ftJ 



(9) 



where h and II are a typical length and a typical 
velocity in the flow and K is a dimensionless parameter 
that is a measure of the ratio of the times in which 
temperature cha:,res occur in the gas to the relaxation 



time of the 



gas 



It will be seen later that is de- 



fined to make it become unity afte^ an instantaneous ex- 
pansion which starts from rest with equilibria energy 
partition and ends with the velocity U. Eliminating T 
between equations (7) and (8) and introducing the non- 
dimensional quantities gives 



ax 



du*' 

at 1 



(io) 



If u T (t T ) is known, the integral of equation (10) can 
be written as 
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i = e 



•/Kdt»/ 



' n du' 2 /Kdt» 
e 



N 



dt 



dt 1 + Constant 



(id 



The rate of entropy increase in the flow can now be calcu- 
lated from equation (11). The rate of heat flow from the 
temperature T vlb to T is k:j hence, 



dS 
dt 



— i 



1 \ 



T. 



rib/ 



(12) 



Now 

as 



e = c 



vib ( T vib " T ) 



- t\ and equation (12) can be written 

1 \ 



dt 



T + 



—J 

vib/ 



The entropy increase along the streamline in question 
between the starting time t and the time t is 



n' 



AS - / 



to/l — 



dt 



(15) 



c vib/ 



In order to obtain the total entropy increase, equa- 
tion (13) would have to be integrated over all the 
streamlines in the flow with the use of equation (7) . 



Similarity Law for Low-Velocity Flows 

The calculation of energy dissipations can be simpli- 
fied if the restriction to flows involving pressure and 
temperature changes that are small compared with ambient 
pressure and temperature is adopted. The greatest ad- 
vantage of this procedure is that the flow pattern ob- 
tained in an incompressible fluid can be used as an ap- 
proximation. This fact is important because few compres- 
sible fluid flows are known accurately. If this 
restriction is accepted, k and hence K can be assumed 
constant for the flow. Equation (11) then becomes 



e« = e 



-Kt' /^du' 2 Kf 



dt 



dt' + Constant 



(14) 
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Now both arid variations of T are small compared. 

Ovib 

with T and equation (13) becomes 

AS m — / e 2 dt (15) 

z o 

It is now shown that there is a simple relation 
among the dependencies of the energy dissipation in a 
low-velocity flow on the scale of the flow, on the typical 
velocity/ and on the relaxation time of the gas. This 
relation is that the entropy increase, reduced to non- 
dimensional form, depends in geometrically similar flows 
on a single "parameter K. • 



Equation (15) can be rewritten as 

\ 2 v JJ t * 



A3- IL_ 1_ J g ' dt ' 



0 



Introducing the nondimensional entropy increase A3 1 by 
dividing A3 by the entropy increase following an 
"instantaneous" compression gives 

AS' s 2K f e' 2 dt' (16) 

* 0 

# 

Prom equation (14), it is kncwm that € 1 and hence AS f 
depends only on K for similar flows. 
> 

Approximations for Large and Small Values of K 

The integrations of equations (14) and (16) are 
sometimes difficult to perform analytically and laborious 
to evaluate numerically. For the special cases in which 
the relaxation time is cither long or short compared with 
the times in which temperature changes take place in the 
gas, it is possible to use approximations that consider- 
ably reduce the numerical labor. In these cases, it is 
possible to express AS* in terms of integrals in which 
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K does not appear under the integral sign; thus, these 
integrals need be evaluated only once to determine AS » 
for all values of K for which the approximation is 
valid. 

The case of short, relaxation time, when K is large, 
will he treated first. Tn order to avoid confusion, the 
symbol t f a is introduced into equation (14), which 

becomes 



£ , (t , a) „ /*** as^-K^-tv 

z 0 



dt 



(17) 



For a large value of K, most of the contribution to 
this integral comes from values of t 1 so close to t' a 
that the following approximations can be made: 



du' 2 _ (du< 2 \ 
dt' \dt' /a 



t' - t' a = 



u' 



u' 



a 



\dt' 



a 



and the lower limit of the integral in equation (17) can 
be replaced by ±°o . Equation (17) then becomes 



c'U'a) = / 

V±ro 



r 



exp 



du' * \ 
, dt ' k 



- u' 



2- 



du 



,2 



where the sign of the lower limit is opposite that 
'du ,2 \ 



of 



dt' 



/a 



Hence, for K » 1, 



^ a/ " K \ dt > } & 



r 



AS' = S 



2 / / du 



,t2 



K i \dt» 



2 



dt' 



(18) 



The case of long relaxation time, when K is small 
compared with 1, is now considered. In the usual flow 
problem, the gas velocity changes appreciably during a 
certain time interval - say, from 0 to - and then 
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settles to a new steady value. The problem can be 
divided into two parts: 0 < t* < t*x and t ! > t f !• 
If K is snail enough, the change in c f due to the 
5 '-term in equation (10) is small compared with the 

ehange due to the — term and can be neglected in 
the calculation of the entropy increase AS ! ]_ during 
the first interval; thus, 

o o 



where u ! q^ is the velocity squared at t ! = 0. Hence, 



1 X 9 



A3» x = 2K I c^dt* = 2i: - a'n*J dt' 



/ 1 / ° "\2 
K \ ^u f " - u»o j at 1 



In order to compute the value of c* at t ! ]_, the 
total contribution of the € t -terra in equation (10) is 
added to the total change in the square of the velocity 

during the first interval £u f °. us, 

«!-. =: AU' 2 - K 



J 0 J 0 ' 

In the period after t f Tj 



€ i = € ? i e 



and the entropy increase in this second period AS'o is 



f 1 

i 
j 



A3' g = 8IC€ ' ^ ' e v J -' dt' = £ ' - ± 



The total entropy increase in the flow is therefore, for 
K « 1, 
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Calculation of Total-Head Defect in Plow 

about a "Source-Shaped" Impact Tube 

The total-head defect to be anticipated in a com- 
pression at the nose of an impact tube of a special shape 
is calculated to be used in the measurement of the re-^ 
laxation time of gases. The restriction to low veloci- 
ties adopted previously is retained, chiefly to permit 
the use of incompressible-fluid theory and of the simi- 
larity theorem. 

The flow about bodies of revolution in a uniform 
stream is usually calculated by considering the flow 
about sources in the fluid. (Compare reference 6, 
p. 146.) It is possible to find a surface in the flow 
across which no fluid flows. If a solid body shaped 
like this surface is substituted for the sources, no 
alteration outside the surface occurs j the flow about 
the solid body is thus identical with that about the 
sources. The flow about a single source in a uniform 
flow is calculated In reference 6 and the corresponding 
shape is plotted in figure 5. The total-head defect 
to be anticipated for a tube of this shape is calcu- 
lated as follows: 

The velocity along the central streamline is re- 
quired. This velocity is given on page 147 of refer- 
ence 6 and is plotted in figure 5 as 

u(x > - 4 - £) ■ 
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where 

x distance along central streamline from source 

U velocity far from body 

d diameter of "impact tube 

This expression can be converted to the following non- 
dimensional form by using U as the typical velocity 
and d as the .typical dimensions 

u i (x ») = ! m _JL_ (19) 
16x f c 

The next step is to find u T (t T )« The quantity t T can 
be found as a function of u r by integrating 

dt« - m - 1 du ' 



u'U') 8 u , (1 I ut) 3/2 



The choice of the zero of t T is arbitrary, 
venience, if t T = 0 when u T = 0.99, then 



For con- 



COl' 



log (f ♦ 
V 9 1 + a/1 - u7 



2 



1 - u« 



- 20 



(20) 



The next step is to determine e'(t') from, equation (14). 
Then, by use of equation (19), 



flu*' 
dt< 



* 2u< 2 Sgi = -16u' 2 (l - u') 3/2 - 
dx 1 v ' 



(21) 



Because e f is zero initially (t* ss -co) and remains zero 
until u f bepins to vary rapidly with tir.e, if K is not 
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too small, the lagging heat capacity can be assumed to 
follow the temperature changes in the gas up to the point 
u T = 0.99; that is, = 0 can be used for t T = 0. 

Combining this fact with equations (14) and (21) yields 

1611*2(1 - u') 5/2 e Kt ' dt< (22) 



In view of the partly transcendental nature of equa- 
tion (20), it was necessary to integrate equation (22) 
numerically. Equation (20) was plotted (fig. 6) in 
such a way*' that the values of u T corresponding to 
regularly spaced values of t T could be found easily. 
By Simpson 1 s rule, e f (t f ) was then found for a series 
of values of K. An example of the result of such a 
calculation is given in figure 6 for K =3. The 
entropy increase along the central streamline was then 
found from equation (16). 

Values of AS * found from integrating equation (22) 
by Simpson 1 s rule and equation (16) with a planimeter 
are plotted in figure 7 and are given in the following 
table : 

RESULTS OP NUMERICAL CALCULATIONS OP AS ! 
FOR SOURCE-SHAPED IMPACT TUBE 



K 


AS' 


10 


0.1685 


3 


.405 


2 


.516 


1 


.676 


.3 


.868 


.1 


.952 



For large and small values of K, the approximations 
developed earlier were used to reduce the labor of cal- 
culations and yielded the result AS' = 1.743/K when 
K is large and AS' ss 1.452K + (1 - 1.008K)* 2 when K 
is small. ' These results are plotted in figure 7; this 
figure thus indicates the range of applicability of these 
approximations . 



e' (f) = -e" 
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Calculation of Entropy Increase in Plow 

through a Nozzle of Special Design 
^ For the measurements of the relaxation time in COo, 

LPs 

~-t a nozzle Is employed in which the gases expand and ac- 

^ celcrate before meeting the impact tube. This expansion 

cannot always be made slow enough - that is, the nozzle 
large enough - that the expansion through the nozzle in- 
volves a really negligible entropy increase; hence, the 
results of figure 7 must be corrected for the entropy 
increases in the nozzle. In order to simplify the cal- 
culations, the nozzle was so designed that the time rate 
of temperature drop was constant. It can be shown that 
the entropy increase in a nozzle of this design is 

where Km- - ~ and TJ is the final velocity attained 
w TIT 

by the flow m the nozzle. It must be remembered that 
the calculations for the impact tube presumed e T to be 
zero initially. This condition is the case only if 
Km » 1 and hence 'the calculation given here is valid 
only for this case. 

Prom the definitions of K and Km, it is seen 

that Ktvt = *r& and hence the total entropy increase 
c 

i tuo e n 



can be expressed as a function of K alone for a given 
l/d. This total entropy increase is plotted in fig- 
ure 8 against K for the two values of l/d used in 
these experiments and for l/d . 



MEASUREMENT OP RELAXATION TIME OF 00 o 

The theory v T ill now be applied to the measurement 
of the relaxation time of COo. This v;ork ras 
undertaken both to test the theory and to develop a 
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technique that would supolement the sonic methods pre- 
viously used for measuring relaxation times. The 
method essentially consists in expanding the gas through 
a known pressure ratio in a nozzle and compressing it 
again at the nose of a source-shaped impact tube. The 
resultant total-head loss is divided by the total-head 
loss that would be obtained in a very slow expansion and 
a fast compression (equation ( 4) ) • This nondinensional 
total-head loss is compared with a theoretical result 
such as is shewn in figure 3 and the value of K appro- 
priate- to the flow is found. From this value of K, 
the .relaxation time of the gas can be easily computed if 
the velocity before compression and the diameter of the 
impact tube are know. 

During the compression of the gas, the temperature 
and pressure rise from T^ and to Tg and pg # 

respectively. The relaxation time and the hoat capacity 
of the gas thus change along a streamline. ■ The pro- 
cedure previously outlined then gives an average relaxa- 
tion time for the flow. It is assumed that this average 
relaxation time is the relaxation time appropriate to 
conditions halfway between compressed and expanded con- 
ditions. Because p is_olose to p^ and Tg £ T Q , 
these conditions j> and T can be found from 



< * 2 

and 

- T 0 + ?1 

The errors introduced in this way certainly are no 
greater than those due to the low-velocity assumption 
introduced in the theory upon which figure 8 is based. 

Gas 

The gas used in these experiments was commercial 
"bone-dry" COg. This gas was dried by passing it 
through calcium chloride and. then deJrydrite while it 
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was at a pressure greater than 40 atmospheres. The puri- 
fication procedure was not so thorough as methods used, 
in some previous investigations, and it is to be expected 
that somewhat shorter relaxation times would be obtained. 
The primary object of this work is to establish the self- 
-sj consistency of this test method rather than to obtain an 

accurate relaxation time for pure COg. 

Apparatus 

The apparatus used is essentially the same as. that 
schematized in figure 2. A longitudinal section through 
a chamber of the most recent design is shown in figure 9. 
(The chamber used in the tests discussed in the next sec- 
tion did not incorporate the liner and the gas entered 
from the bottom.) The gas enters through three holes 
that were made small to stabilise the gas flow into the 
chamber. The glass v:ocl is necessary to remove turbu- 
lence from the gas in the chamber and contributes materi- 
ally toward reducing the total-head defects obtained in 
gases without heat-capacity lag. It was found that 
total-head defects traceable to nonunif ormi ties in tem- 
perature existed. and could be reduced by the use of the 
lined chamber shown. The fact that the gas flows around 
the inner chamber before entering helps to keep the gas 
in the inner chamber at uniform temperature. 

The temperature nonunif ormi ties can be almost elimi- 
nated if the gas entering the outer chamber is at the 
same temperature as the chamber* A mechanism was used 
to adjust the alinement of the impact tube without moving 
the tip from the center of the nozzle. The impact tube 
must be adjustable in. order that small errors in shape 
near the hole will not give spurious total-head defects 
(in helium, for example). The gas and the chamber were 
heated electrically and a thermocouple inside the chamber 
was used to measure the gas temperature. 

The nozzle used had a circular cross section, was 
1.6 inches long, and was designed according to the methods 
previously described to give a constant time rate of 

du^ 

temperature drop; that is, - Constant for the first 

dt 

1.5 inches, the last 0.1 inch being straight. The 
radius of the nozzle r Is plotted against the distance 
along the center line x in figure 10. 
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Two impact tubes with diameters 0.0299 inch and 
0.0177 inch were used in these experiments. They were 
made by drawing out gifts s tubing until a piece of ap- 
propriate diameter and hole was obtained. The hole was 
kept larger than about 0.004 inch and the fine section 
not too long (~l/4 in.) to prevent the response of the 
alcohol manometer from being too sluggish. The ends 
of the tubes were ground to a source shape (fig. 5) on 
a fine stone. During the grinding process, a silhouette 
of the tube was cast on the" screen of a projecting micro- 
scope and the contour superimposed on a source-shaped 
curve. By this technique the contour could be ground 
to the source shape within 0.0005 inch, except for the 
hole, in a short time. 



Tests and Computations . 

The total-head defect in COg was measured with the 
two impact tubes over a range of chamber pressures. The 
consistency of relaxation times obtained at various pres- 
sure ratios and with various impact tubes serves as a 
check on this method of measuring relaxation time and on 
the theory on which the method is based. 

Before each series of measurements nitrogen, which 
has only a negligible vibrational heat content at room 
temperature, was run through the chamber to be sure that 
no spurious effects and leaks were present. In the re- 
sults reported herein, the errors due to these effects 
were kept to less than 0;01 percent of the chamber pres- 
sure: therefore, the resultant error in * relaxation time 
due to these causes was less than 4 percent. In sub- 
sequent work (not reported herein) , it was found that 
most of these total-head aberrations could be eliminated 
by ensuring uniform temperature - in the issuing gases. 
If care is taken to eliminate temperature nonunif ormities , 
tube misalinements, and turbulence in the chamber,* the 
total-head aberrations can be reduced to 0.002 percent or 
less . 

The total-head defects obtained were divided by the 
result of equation (4) to reduce them to nondimensional 
form. The appropriate value of K was found by refer- 
ring to the appropriate curve in figure 8. The gas 
velocity was computed from the reading of the mercury 
manometer by the enthalpy theorem with adiabatic expansion 
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assumed. The relaxation tine was then computed from the 
definition of K by equation (9) . The relaxation times 
thus obtained were expressed in collisions per molecule. 
The number of molecular collisions per second in COg 

was assumed to be- 8-888 x 10 9 at 15° G and 1 atmosphere 
by combining tables of pages 26 and 149 of reference 7. 
At all other temperatures and pressures, the number of 
collisions was assumed to vary inversely with y/T and 
directly with pressure. The number of^molecular colli- 
sions per second and the heat capacity of the gas were 
computed at temperature T and pressure p. The data 
obtained are given in tables I and II for the 0.0299- 
and 0.0177-inch tubes, respectively. 

The results are plotted in figure 11, which indi- 
cates that the relaxation time in collisions is nearly 
independent cf pressure ratio and impact-tube size. 
This consistency constitutes the desired verification 
of this test method. It was expected that a variation 
at high pressure ratios would appear in view of the 
assumption of low velocity made at several points in the 
theoretical development . 

A large part of the scatter of the results in fig- 
ure 11, in particular the apparent drop at low pressures, 
is attributed to the fact that in the tests the average 
temperature (halfway between chamber and expanded tem- 
peratures) was not held constant during the run. 

The average number of collisions obtained with the 
0. 0299-inch tube was 33,100; with the C. 0177-Inch tube, 
32,000. The final result at 105° ? thus is 32,600, 
which is somewhat lower than the result of recent inves- 
tigations in which the COg has been much more carefully 

purified than in the present investigation. (Compare 
with fig. 1.) 



IMPACT -TUBE METHOD OP MEASURING RELAXATION TITS OP GASES 



The impact-tube method of measuring the relaxation 
time of gases rests essentially on the fact that the 
total-head defect not traceable to heat-capacity lag can 
be reduced to a very small value - say, 0.002 percent. 
Very small dissipations due to heat-capacity lag are 
therefore measurable. For example, a gas having a 
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lagging heat capacity 0.1R with a relaxation time of 10""' 
second could give a total-head defect of 0,05 percent. 

If the gas had a lagging heat capacity as large as R, a 

p 

relaxation time as short as 10""° second would be meas- 
urable • 

This method seems to be easier to carry out than the 
sonic methods previously discussed and can be used to 
measure relaxation times with comparable precision. The 
quantity of gas required to make a measurement will be 
larger than for the sonic methods (a standard tank of CO2 
lasts about 5 hr in this apparatus) and thus may make it 
more difficult to attain high purity. 

If the gUS to be studied has a long relaxation time - 
greater than 50 microseconds, for example - it should be 
possible to measure the relaxation time in an apparatus 
similar to the one discussed by comparing the total-head 
defects obtained with a calculation of the entropy in- 
crease in the nozzle. In this case, the time taken for 
the gas to flow through the nozzle is compared with the 
relaxation time of the gas. The shape of the impact 
tube would be unimportant in this case as long as it was 
small enough that K « 1. 

CONCLUSIONS 



The existence of energy dissipations in gas dynamics, 
which must be attributed to a lag in the vibrational heat 
capacity of the gas, has been established both theoreti- 
cally and experimentally. 

An approximate method of calculating the entropy in- 
crease in a general flow problem has been developed. The 
special case in which a gas at rest expands out of a 
specially shaped nozzle and is compressed at the nose of 
a source-shaped impact tube near the mouth of the nozzle 
has been treated, and the dependence of the resultant 
total-head defect on the relaxation time of the gas has 
been found. 

The total-head defect in this flow has been applied 
to measure the relaxation time of C0 2 . The results 
obtained with two impact tubes were in agreement within 
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about 3 percent. The consistency of these results is 
regarded as a check on the general theory developed and 
on this measurement method. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. 
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APPENDIX A 
SONIC MEASUREMENTS IN CARBON DIOXIDE 

Much careful work has been done on the lag in the 
vibrational heat capacity of CCg. Carbon dioxide is a 
linear molecule and thus has a transnational and rota- 
ticnal heat capacity of ~R. It has four normal modes 

in vibration that are diagrammed with their frequencies 
a3 follows (data from reference 8) 2 

0 — > C <— 0 y) ± » 4,164 x 10 13 
0 C 0 u 2 = 2c 003 x 10 15 2 modes 



I 4- 



C — > < — 0 v~ - 7*050 x 10 



13 



The heat capacity of COg is somewhat complicated 

by the fact that the second excited state of the oscil- 
lation u 2 has almost the same energy as the first 
excited state of ^ . The near resonance results in a 
strong interaction through the first-order perturbation 
(the first-order departure of the potential energy from 
the square law) between the two states involved, as was 
pointed out by Fermi (reference 9). This perturbation 
produces significant disturbances (*5Q cm*- 1 ) of the levels 
involved but does not have a large effect on the heat 
capacity of the gas. The heat capacity of CO2 was 
computed by Kassel (reference 10) and his results are 
used in the present calculations . 

Eucken and his coworkers have carefully studied over 
a period of years the dispersion of sound in COg 
(references 11 to 15). One conclusion of this work - 
that the vibrational energy levels in COg adjust with 

the same relaxation time - is demonstrated by showing that 
the dispersion curves obtained fit a simple dispersion 
formula such as Kneser ! s . 

Ku.fcfcLLfccr:, for example, obtained a simple dispersion 
curve at 410® C, at which appreciable heat capacity due 
to all three normal modes would be expected. Richards 
and Reid (reference 16) and others (see bibliography of 
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reference 2) have maintained that the symmetrical valence 
vibration u-i of C0 o does not adjust at 9 kilocycles 
in some dispersion measurements made near room tempera- 
ture. As they point out, this fact Is remarkable be- 
^- cause the second excited state of u g strongly perturbs 

the first excited state of u n . In any case, the con- 

tribution of this normal mode to the heat capacity is 
very mall at room temperature and the effects found are 
near the limits of the accuracy of Richards and Reid* 

The relaxation time of C0 2 in molecular collisions, 
as ^given by Eucken and his coworkers, is plotted against 

T 3 (T in °K) in figure 1 for comparison with the theory 
discussed in appendix 13 ♦ Van Itterbeek, de Bruyn, and 
Mariana (reference 17) measured the absorption at 599 kilo- 
cycles in very carefully purified COg. Their measure- 
ments, which are also given in figure 1, show a longer re- 
laxation time than the measurements of Eucken and his 
coworkers. They attribute this increased relaxation time 
to careful purification of the gas. 

All the measurements with C0 2 have indicated that 
the relaxation time is inversely proportional to the 
pressure of the gas. This result shows that the process 
responsible for the interchange of energy between vibra- 
tional and other decrees of freedom is bimolecular. 



APPENDIX B 

THEORY OF EXCITATION OP MQUHWIAH VIBRATION BY COLLISIONS 

Lardau and Teller (reference 18) have given an ap- 
proximate calculation of the probability of the excita- 
tion of a- vibrational quantum in a molecular collision. 
Assuming that the interaction energy which induces the 
vibration depends linearly on the nam*! coordinate of 
a harmonic vibration, they make a first-order perturba- 
tion calculation. The matrix element for the transi- 
tion from the 2 th to the (I + l)th or from the 
(I f 1) th to the Ith vibrational state is then pro- 
portional to \A " +"i . The transition probabilities 
k i.i are proportional to the square of matrix elements 
and therefore 

k 2 k ik 53 V fir o lr — 1 ■ o • ^ 

01* 12- 23 ^10 ok 21- K 32 l.rfto 

and, when i - j ^ +1, i z± . - 0. This result is shov;n 

in reference 18 to lead to the prediction that all the 
allowed transitions in a given normal, mode have the same 
relaxation time. Landau and Teller next examine the 
collision process classically, assuming the interaction 
energy between translation and vibration to be propor- 
tional to 0 r 9 where x is the &t stance between the 
molecules and a is an undetermined constant. They 
also assume that the trans la tional energy of the mole- 
cules - that is, the collision - is adiabatic. The 
amount of energy transferred to vibration in a collision 
is then calculated and used to estimate the transition 
probability k 01 and the relaxation time of the gas. 

They conclude that the temperature variation of the re- 
laxation time expressed in molecular collisions is given 
by 



Collisions 




where M is molecular weight • 

In figure 1, experimental results for the relaxation 
time In collisions of C0 2 and nitrous oxide NgO are 
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JL 

plotted against T 3. The theoretical results are seen 
to be straight lines, within experimental error. The 
value of a can be found from the slope of the straight 

line. For CO2 with u = 2.003 x 10 13 , a = 0.22xlO~ S cm 

and, for KfO with 0 b 1.773 x 10 15 , a=0.36xl0~ 8 cm. 

These reasonable values for a are a further check on 
this theory. 

It should be pointed out that the temperature varia- 
tion of catalytic effects is quite different from that of 
pure gases, the number of collisions required being nearly 
independent of temperature. (See Kuchler, reference 15.) 
Various attempts have been made to associate the effective- 
ness of catalysts with their physical or chemical proper- 
ties but no generally successful rule seems to have been 
proposed. Gases that have some chemical affinity, gases 
with large dipole moments, and gases with small moments 
of inertia are usually most effective. 
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APPENDIX C 
THE ENTHALPY THEOREM 

If no energy is transmitted across the walls of a 
stream tube, the total energy (internal energy E plus 

i ah 

Kinetic energy per unit mass *- u / plus the work done by 

pressures pV must be the same at any cross section of 
the tube ; that . is, 

% ' -j p 

E + ~u + pV = Constant (CI) 

o 
m 

In the case of a perfect gas with constant heat capacity 
and with equilibrium partition of energy, equation (01) 
becomes 

i 2 

c T + -u w = Constant (02) 
P 2 

where c p is the heat capacity at constant pressure and 

T is the absolute temperature* Whenever equilibrium 
partition exists, even though nonequilibrium states have 
been passed through, equation (01) is applicable in the 
absence of viscosity and heat conduction and equation (02) 
can be applied to perfect gases, " provided the heat 
capacity" of the gas can be considered constant. 
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TABLE I.- DATA AND CALCULATIONS FOR 0.0299-INCH TUBE 



2 



8 



Chamber 
pressure, 
' P 0 / Pl 

(atm) 



\otal-head 
defect, 



£o 

Pp. 



- 1 



100 

p 0 

— - 1 

(percent) 



Average 
tempera- 
ture, 
f 

(°P) 



Total-head 

defect 
from inst. 
compression 
theory, 

PO 



- 1 



Total 

A3' 
/ Col. 2 \ 
V'Col. 4/ 



I 

(from 
fig. 3) 



loo- 



t'2 



Gas 
velocity 
at end of 
nozzle , 
U 

( in. /sec) 



~- - 1 

Pi 

(percent) 



Relaxa- 
tion 
time, 

r= JL 

Uk 

(r.iicro- 
sec ) 



Collisions 
(Y x Collisions 
per molecule 
per sec) 



1.185 
1.253 
1.323 
1.388 
1.464 
1.527 
1.580 
1.672 
1.734 
1.738 
1.295 
1.350 
1.364 
1.422 
1.499 
1.552 



0.321 
.482 
.669 
.805 
.966 
1.088 
1.193 
1.344 
1.445 
1.467 
.586 
.721 
.758 
.880 
1.025 
1.130 



115.0 
109.8 
107.5 
ICO. 7 
101.6 
98.2 
96.1 
93.9 
92.6 
92.4 
107.8 
105.4 
111.5 
107.2 
99.2 
98.1 



0.611 
.815 
1.024 
1.198 
1.412 
1 . 577 
1.706 
1.934 
2.076 
2.084 
.941 
1.098 
1.152 
1.312 
1.502 
1.641 



0.525 
.591 
.653 
.672 
.684 
.690 
.699 
.695 
.696 
.704 
.623 
.656 
. 658 
.671 
.682 
.688 



2.03 

1.58 

1.23 

1.13 

1.07 

1.05 

1.01 

1.025 

1.02 

.982 
1.38 
1.21 
1.21 
1.14 
1.075 
1.05 



5640 
6464 
7188 
7731 
8336 
8761 
9090 
9611 
9936 
9955 
6910 
7423 
7591 
8051 
8573 
8926 



2.61 
2.93 
3.38 
3.42 
3 • 35 
o • 2 5 
3.26 
3.04 
2.95 
3.06 
4.17 
3 • 33 
3.25 
3.26 
3.25 
3.19 



JO 

35 



24,000 
27,800 
33,200 
34,700 
35,100 
, ,000 
35,900 
34,700 
34,500 
35,900 
30,400 
33,100 
32,400 
33,300 
34,500 
34,700 



w. 33,100 



TA3LTH II.- DATA AND CALCULATIONS FOR 0. 0177-INCH TUBE £ 



1 


2 


3 


4 


5 


6 


7 


8 


9 


Chamber 
pressure , 

P 0 /Pi 
(atn) 


Total-head, 
defect, 

P 0 
P2 

100 

Po 1 

Pi ' 

(percent) 


Average 
LGrnpoi a— 
ture, 

(°P) 


Total-head 

defect 
from inst. 

c nnuopfis s ion 
theory, 

:°2 

100- 

^ - 1 

Pi 

(percent) 


Total 
A3' 
(Col. 2\ 
VCol. 4/ 


K 

( from 
fig. 8) 


Gas 
velocity 
at end of 
nozzle , 
U 

( in. /sec) 


Relaxa- 
tion 
time , 

d. 

T= Uk 
(micro- 
sec) 


Collisions 
(T x Collisions 
per molecule 
per sec) 


1.442 
1.315 
1.334 
1.389 
1.278 
1.229 
1.205 
1.467 
1.686 
1.612 
1.665 

JL • vJ CJ JL 

1.556 
1.548 
1.531 
1.515 
tL.482 
1.417 


1.014 
.765 
.748 
.900 
.633 
.519 
.430 
1 . 100 
1.521 
1.383 
1;495 
1.317 
1.283 
1.276 
1.244 
1.216 
1.123 
.986 


97.7 
127.3 
106.1 
98.8 
100.8 
110.1 
121.1 
103.4 
95.2 
96.0 
95.0 
95.1 
96.1 
99.2 
100 . 9 
104.3 
102.9 
107.3 


1.342 
1.041 
1.057 
1.196 
.876 
.742 
.683 
1.429 
1.972 
1.788 
1.919 
1.708 
1.644 
1.637 
1.594 
1.564 
1.468 
1.298 


0.756 
.735 
.708 
.753 

.699 
.704 
.770 

« ( f C* 

.774 
.779 
.771 
.730 
.780 
.781 
.778 
.765 
.760 


0.73 
.81 
• 9^ 
.74 
.80 
.97 
.93 
.67 
.67 
.67 
.65 
.67 
.64 
.64 
.64 
.65 
.69 
.71 


8146 

7232 

7289 

7723 

6687 

5182 

5940 

8377 

9703 

9283 

9585 

9088 

8937 

8912 

8807 

8728 

8482 

8013 


2.98 
3.02 
2.64 
3.11 
3.11 
2.95 
3.20 
3.15 
2.72 
2.85 
2.84 

3 . 09 
3.10 
3.14 
3.12 
3.00 
3.11 


31,700 
29,700 
26,700 
32,400 
30,800 
28,400 
30,200 
33,700 
32,000 

33,100 
32,800 

34,400 
34,600 
34,000 
32,600 
32,500 

Av. 32,000 
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Figure 1 Measurements of relaxation time of CO2 and N2O. 
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Figure 4.- Instantaneous compression, theoretical and experimental. 
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Figure 7 The nondimensional entropy increase in the flow approaching a source-shaped 

impact tube*. " K = d/UT. 
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Figure 9 • - Longitudinal section through chamber. 
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Figure 10.- Shape of no: "zle t radius against dis- 
tance along center line. 
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Figaro 11..- Ivlea surer:: en t of relaxation ti;:ie of vib: 
105° ?. Collisions = Relaxation time > 



f vibrational heat capacity in CO^. Temperature near 
x Collisions per r.olecule per second. ' 
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